WHY SHOULD ONE EXPECT TO FIND LONG RUNS OF 
(NON)-RAMANUJAN PRIMES ? 



PETER HEGARTY 

ABSTRACT. Sondow et al have studied Ramanujan primes (RPs) and observed nu- 
merically that, while half of all primes are RPs asymptotically, one obtains runs of 
consecutives RPs (resp. non-RPs) which are statistically significantly longer than one 
would expect if one was tossing an unbiased coin. In this discussion paper we attempt 
a heuristic explanation of this phenomenon. Our heuristic follows naturally from the 
Prime Number Theorem, but seems to be only partly satisfactory. It motivates why one 
should obtain long runs of both RPs and non-RPs, and also longer runs of non-RPs 
than of RPs. However, it also suggests that one should obtain longer runs of RPs than 
have so far been observed in the data, and this issue remains puzzling. 



1. The model 

Consider the following random process : you have an infinite supply of identical 
biased coins, which return heads with probability p G (|, 1], and tails with probability 
1 — p. Now toss iV of these coins. For each i = 1, N, let hi, ti denote the number of 
heads (resp. tails) among the first i tosses. Thus hi + ti — i. Also denote 

Ai:=hi-ti. (1.1) 

The coins that come up as heads will each be colored red or blue, according to the 
following rule : Suppose the i:th toss is a head. Then color this coin red if and only if 
the following two conditions sare satisfied : 

Aj > A i; fora\\j=i + l,...,N (1.2) 

and 

If 1 < k < i and A/. > Aj, then there exists I E (k, i) such that A^ < Aj. (1.3) 

Otherwise, color the coin blue. The definition requires some thought, so here is an 
example to illustrate how the scheme works : 
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Now what do we expect to observe, when N is large ? Well, with high probabil- 
ity (w.h.p.), we will observe close to pN heads and close to (1 — p)N tails. Hence 
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Ajv ~ {2p — l)N w.h.p. Now at most one coin is colored red for each positive value 
attained by the function A. On the other hand, for any e > 0, w.h.p. the function A will 
eventually exceed (2p — 1 — e) N for good. Hence, w.h.p. about {2p — 1)N coins will be 
colored red, and the remaining heads, about (1 — p)N in number, will be colored blue. 
So, as iV — > oo, the fraction of redheads, amongst all heads, will almost surely (a.s.) 
approach and the fraction of blueheads will a.s. approach In particular, when 
p = 2/3, about half the head-coins will be colored red and half colored blue. 

I claim that, with p = 2/3, this is a good basic model to have in mind when one 
considers Ramanujan primes (RPs) : the red coins corresponding to RPs and the blue 
ones to non-RPs. I will explain two things : 

1. Why this model is reasonable. 

2. Why one expects to get longer runs of blue coins than if the red-blue coloring was 
done by tossing another, fair coin. 

I will deal with the second issue first. However, the analysis will show that, in this 
model, one also expects longer runs of red coins than if the coloring was done fairly at 
random, though not as long as the blue runs. This may seem to contradict the data in 
HSNNH . After explaining why I nevertheless consider the model to be reasonable, I will 
discuss this issue. 



2. Why do we get long monochromatic runs ? 

If a biased coin with probability p of heads is tossed N times, then it is well-known 
that the expected length of the longest run of consecutive heads is approximately iog(i/p) = 
log 1 / p N. Another way of looking at this is that, for any k G N, one expects to have 
to toss the coin on the order of (l/p) k times to have a reasonable probability of seeing 
at least one run of k consecutive heads. This is easy to see intuitively : the probability 
of any k consecutive coin tosses all resulting in heads is p k and thus, by linearity of 
expectation, the expected number of such runs amongst N tosses is (N — k + l)p k , 
which (for any fixed k) will be 0(1) when iV = Q[(l/p) k ]. 

In particular, when p = 1/2, we expect to have to make on the order of 2 k tosses to 
have a reasonable probability of witnessing a run of k heads. 

The following facts about biased coin-tossing are also well-known : 

Proposition 2.1 Suppose we toss a sequence of identical biased coins with probabil- 
ity p > 1/2 of heads. With notation as in Section 1, for each N e N, let cn, p denote the 
probability that Aj > Ofor alii = 1, N. Then the numbers cn, p are non-increasing 
in N and if we let c p := limAr^oo cn, p , then 

2p- 1 

c«=— >0. (2.1) 

p 

PROOF : That cn. p > cn+i :P is trivial. Let E, F and F' denote the following three 
events : 
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E : the event that A» > for all i > 0, when we make an infinite sequence of tosses. 
F : the event that Aj > for alH > when we make an infinite sequence of tosses. 
F' : the event that Aj > Ai for all i > 1 when we make an infinite sequence of tosses. 

Since the coin tosses are independent, one has 

P(F) = P(F'). (2.2) 

Secondly, it is clear that 

Cp = P(F). (2.3) 
Thirdly, the event E occurs if and only if the first toss yields a head and thereafter event 
F' occurs. Hence, 

P(E)=p-F(F'). (2.4) 
In Example 1.13, Chapter 3 of (D), it is shown that 

F(E) = 2p-l. (2.5) 
Eqs. (2.2)-(2.5) together imply (2.1), and the proof is complete. 

Remark 2.2 It is no accident that c p equals the fraction of redheads in the head-coloring 
model of Section 1. Indeed, this observation is the basis for the rigorous proof of (2.1). 

Fix p > 1/2 and consider the head-coloring model of Section 1 . Fix k G N and let 
Et,N denote the expected number of runs of k redheads, when N coins are tossed. I 
claim that, as iV — > oo, 

P*" 1 < %^ < P*- 1 . (2-6) 
p J JM 

To see the right-hand inequality, just observe that if k consecutive heads are all colored 
red, then at the very least the k — 1 heads from the 2nd to the last must have been a run 
of k — 1 heads, with no tails in between. This happens with probability p k ~ l . Thus, the 
probability of a run of k redheads with a fixed starting point is at most p k ~ l . Since there 
are N — k + 1 possible starting points for the run, linearity of expectation implies that 

£ k , N <(N-k + l)p k -\ (2.7) 

which gives the right-hand inequality in (2.6). For the lower bound, we again consider 
a fixed starting point. A sufficient condition to get a run of k redheads with a given 
starting point is that the following three events all occur : 
A : the starting point is a redhead, 

B : it is followed by a run of k — 1 heads, with no tails in between 

C : the value of A never again goes below its value at the end of this run of k heads. 

It is clear that each of B and C positively correlates with A, while B and C are in- 
dependent of one another. Hence 

¥(AAB AC) > F(A) x ¥(B) x P(C). (2.8) 

As shown in Section 1, we know that ¥(A) — > 2p — 1 as iV — > oo. As above, the 
event B occurs with probability p k ~ l . Thirdly, it is immediate that P(C) > c p = 22^. 
Plugging everything into (2.8), we find that the probability of a run of k redheads with 
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a given starting point is at least ( p ' -p k l . Linearity of expectation then yields the 

left-hand inequality in (2.6). 

This brings us to our first result : 

Proposition 2.3 In the model of Section 1, the expected length of the longest run of 
consecutive reds among a total of N heads, is on the order of ^^n^\ ■ In other words, 

for any k G N, we expect to have to make on the order of p k tosses in order to have a 
reasonable probability of observing a run ofk redheads. 

In particular, when p = 2/3, the expected length of the longest run of consecutive 
reds among a total of N heads, is on the order of 1^(3/2) • ^ n otner wor ds, for any 

k E N, we expect to have to make on the order of (|) k tosses in order to have a reason- 
able probability of observing a run ofk redheads. 

Remark 2.4 The proposition says that, for any fixed p > 1/2 and very large N, we 
expect to see runs of redheads amongst the heads of similar length to runs of heads 
amongst all the coins. 

So what about blues ? Here, for simplicity, I only consider the case p = 2/3 for 
the momenf]. Let k 6 N and let J^2k,N denote the expected number of runs of 2k con- 
secutive blue heads somewhere amongst the first N coins. Suppose, for example, that 
in a run of 3k consecutive tosses one observes at least 2k tailsQ This means that the 
function A will have decreased by at least k over this run. All succeeding heads will 
definitely be colored blue at least until the function A has risen by k again. Then it is a 
tedious, but standard, calculation to show that there is a fixed u > C0 such that, in order 
for the function A to increase by k, at least 2k heads will need to be revealed. Hence, 
as iV — > 00, 

^ > u ■ g 2k , (2.9) 

where q 2 k is the probability of a run of 3k tosses yielding at least 2k tails. Explicitly, 
one has 

3k \ I l \ 1 2 



l=2k 



A lower bound for this is got by simply taking the I = 2k term, hence 

£)(§)"(§ 



'i will generalise to arbitrary p > 1/2 when I get the time. I will indicate below where changes need 
to be made. 

More generally, one will need to replace 2 and 3 by some numbers depending on p, and chosen in 
such a way that the final exponent in (2. 13) will be less than p. 
This number will also depend on p in a general analysis. 
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Now let k — > oo also. Put v k := \f Applying Stirling's estimate to (2.11), one 
easily computes that 



A- 



and hence that 

q2k>v k \^J =v k \^=J . (2.13) 

Putting all this together, and using the fact that u in (2.9) is a constant, plus that the 
function v k decreases subexponentially in k, we have our second result : 

Proposition 2.5 In the model of Section 1, with p = 2/3, the expected length of the 
longest run of consecutive blues among a total ofN heads, is at least lc |°|y^ = 2 log 2 N. 
In other words, for large but fixed k E N, we expect to have to make at most on the order 
of(\^2) k = 2 fc / 2 tosses in order to have a reasonable probability of observing a run of 
k blueheads. 

Remark Note the use of the words 'at least' and 'at most' in Proposition 2.5, as against 
'approximately' in Proposition 2.3. This reflects the fact that we only have a lower 
bound in (2.9), whereas in (2.6) we have both upper and lower bounds. While it seems 
difficult to compute exactly, it is quite easy to see, with the help of Proposition 
2.1, that there will be some upper bound of the form cf k , for some C\ < 1. Hence, 
in order to have a reasonable probability of observing a run of k blueheads, one does 
expect to have to make a number of tosses which is exponential in k. 

From Propositions 2.3 and 2.5 it follows that one expects to see longer runs, both of 
reds and blues, than if the coloring was done fairly at random, but that one expects to 
see even longer runs of blueheads than of redheads. Or, to put it another way, for any 
fixed, and large enough k E N, one expects to see a run of k blueheads somewhat earlier 
than a run of k redheads, and one expects to see both in turn much earlier than if the 
coloring was done fairly at random. 

3. Why is this a reasonable model for RPs ? 
First of all, we just focus on the 'ordinary' RPs discussed in [SNN 1 which we wish 



to compare with the p = 2/3 model above. The obvious extension to the so-called 
Generalised RPs of [ABMRS | will be mentioned at the end. 

So consider p = 2/3 as fixed for now. Let me describe, in somewhat informal terms, 
another random model which I claim is asymptotically equivalent to that in Section 1. 
Suppose we have two different radioactive substances, H (head) and T (tail). H decays 
twice as fast as T (i.e.: T has double the half-life of H). At some time t = 0, I start 
observing both substances and record each individual decay of a H- or a T-atom. For 
each i E N, let hi, U denote the number of H- (resp. T-) atoms amongst the first i which 



4 Table 1 of this paper contains some errors, but these have been corrected in arXiv: 1 105.2249'v2) 
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decay. Then color a decayed H-atom red if conditions (1.2) and (1.3) hold, otherwise 
blue. 

If one wants to be more formal, one can phrase this in terms of two independent, 
parallel Poisson processes, one of which has double the intensity of the other. But the 
point is that this model is equivalent, as t — > oo, with that of Section 1. I leave it to the 
reader to convince himself of this. 

From here, we can see the relevance to Ramanujan primes. First of all, the prime 
number theorem says that, for large x, there are approximately primes up to x. 
Equivalently, it says that, for large n, the n:th prime satsifies p n ~ nlogn. Now sup- 
pose that we start from some p n and begin searching to the right of it for the next 
prime p n +i- There is, of course, nothing random about this search process. But in a 
well-known 'random model' for prime gaps, one imagines that this search is a Poisson 
process with intensity r^— ~ — . In this random model, there is also no need to start 

tr J logn logp„ ' 

the search at a prime : the starting point can also be chosen at random without affecting 
the model. 

Now suppose one chooses a large random number x and starts two prime searches in 
parallel, one at x and the other atx/2, where the former search proceeds twice as fast 
as the latter. By this I mean that, when in the former process we have searched from 
x up to x + t, then in the latter we have searched from x/2 up to (x + t)/2. Since 
log(x/2) ~ log x, one sees that the former Poisson process has approximately twice the 
intensity of the latter. Clearly, the method of determining whether a prime is Ramanu- 
jan or not is basically the same as that of deciding how to color the primes 'revealed' 
in the former of these two Poisson processes. Hence, we have shown how our model in 
Section 1 is a reasonable model for the Ramanujan primes. 

Finally, we turn to Generalised RPs. To model c-Ramanujan primes, one should choose 



Everything above is, of course, heuristics. There is nothing 'random' about the prime 
numbers. More importantly, the Prime Number Theorem is a very precise statement 
about the density of the primes. This suggests that a main problem with our heuristic 
model is its Markovian nature. In other words, if I start a search for a prime from 
some point x, and don't find any prime up to x + t, say, then the PNT implies that I 
am now, in some sense, 'more likely' to find a prime between x + t and x + 2t. The 
further one searches the more restrictive the PNT becomes. Since, in order to observe a 
monochromatic run (equivalently, a run of (non)-RPs) of length k, one expects to have 
to search in a range exponential in k (Props. 2.3 and 2.5), the PNT will, indeed, impose 
severe restrictions on what can happen. 

In HSNNI . the authors fix an upper bound x, and find the length of the longest run 
of (non)-RPs in [l,x]. Our models suggest that a better way of collecting the data 
would be to fix an integer k, and determine the first time a run of k (non)-RPs appears. 
This is because, 'locally', the model of two parallel Poisson processes seems to be okay. 

5 I will leave all further calculations for general c to another forum. 




4. Discussion 
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Despite the problems discussed above, I do not see clearly why Propositions 2. 1 and 2.2 
should not be a reasonable guide what should be observed, as k — > oo. In other words, 
I conjecture that, once the numbers get big enough, one will also observe considerably 
longer runs of Ramanujan primes than a fair coin-tossing model would suggest, though 
never as long as the runs of non-Ramanujan primes. If this is not the case, if the data 
in HSNNB is a reasonable guide to what happens asymptotically, then our model must 
contain a serious flaw. 
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